Spectroscopy of a driven solid-state qubit coupled to a structured environment 
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We study the asymptotic dynamics of a driven spin-boson system where the environment is formed 
by a broadened localized mode. Upon exploiting an exact mapping, an equivalent formulation of 
the problem in terms of a quantum two-state system (qubit) coupled to a harmonic oscillator which 
is itself Ohmically damped, is found. We calculate the asymptotic population difference of the 
two states in two complementary parameter regimes. For weak damping and low temperature, a 
perturbative Floquet-Born-Markovian master equation for the qubit-oscillator system can be solved. 
We find multi-photon resonances corresponding to transitions in the coupled quantum system and 
calculate their line-shape analytically. In the complementary parameter regime of strong damping 
and/or high temperatures, non-perturbative real-time path integral techniques yield analytic results 
for the resonance line shape. In both regimes, we find very good agreement with exact results 
obtained from a numerical real-time path-integral approach. Finally, we show for the case of strong 
detuning between qubit and oscillator that the width of the n-photon resonance scales with the n-th 
Bessel function of the driving strength in the weak-damping regime. 

PACS numbers: 03.65.Yz, 03.67.Lx, 74.50.+r, 42.50.Hz 
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I. INTRODUCTION 



Currently, we witness an impressive progress in real- 
izing coherent quantum dynamics of macroscopic solid 
state device^iiSiSi^. Very recently, experimental results 
on the quantum dynamics of a superconducting flux 
qubit coupled to a read-out SQUID have been reported^. 
The flux qubit consists of a superconducting ring with 
three Josephson junctions and, in the proper parameter 
regime, it forms a quantum mechanical macroscopic two- 
state system (TSS). An external time-dependent driving 
force controls the state of the TSS. A SQUID couples 
inductively to the qubit and, together with an exter- 
nal shunt capacitance, it can be modeled as a harmonic 
oscillator (HO). Due to the coupHng of the SQUID to 
the surrounding environment, the harmonic oscillator is 
(weakly) damped. The state of the qubit can be inferred 
from the state of the SQUID. The experiment provides 
spectroscopic data on the different transition frequencies 
of the coupled TSS-HO device. Moreover, Rabi oscilla- 
tions involving different pairs of quantum states of the 
device have been revealed, including the so-termed red 
and blue sideband transitions between energy states of 
the coupled TSS-HO system. 

Heading for a comprehensive detailed understand- 
ing, a quantitative modeling which includes the effects 
of time-dependent driving, decoherence and dissipation 
is required. Our description goes beyond the well- 
known Jaynes-Cummings model^, by avoiding the strong 
rotating-wave approximation and by including a micro- 
scopic model for the environment. A generic theoretical 
model for studying the environmental effects on a driven 
TSS is the driven spin-boson model^'® where the TSS 
tunneling splitting is denoted by A. The environment 
is characterized by a spectral density J{ijj). The widest 



used form is that of an Ohmic spectral density, where 
J(io) is proportional to the frequency lo. It mimics the 
effects of an unstructured Ohmic electromagnetic envi- 
ronment. In the classical limit this leads to white noise 
and all transitions in the system are damped equally. 
However, if the environment for the qubit is formed by 
a quantum measuring device which itself is damped by 
Ohmic fluctuations, the simple description as an Ohmic 
environment might become inappropriate. In particu- 
lar, the SQUID-detector being well described as a HO 
can equally well be considered as a (broadened) local- 
ized mode of the environment influencing the qubit as 
the central quantum system. In this picture, the plasma 
resonance at frequency Q.p of the SQUID gives rise to a 
non-Ohmic effective spectral density Jeff(a;) for the qubit 
with a Lorentzian peak at the plasma frequency of the 
detector^'^. 

The effects of such a structured spectral density on 
decoherence have been investigated in several theoret- 
ical works in various limits. The role of the external 
driving being in resonance with the symmetric TSS at 
zero temperature has been studied in Ref. ^3 within a 
Bloch-Redfield formalism being equivalent to a pertur- 
bative approach in J^s. Smirnov's analysisi^ is based on 
the assumption of weak interaction between the TSS and 
the HO being equivalent to a perturbative approach in 
Joff as well. Moreover, a rotating-wave approximation 
is used. The first assumption, however, might become 
problematic if the driven TSS is i n reson ance with the 
HO. The results presented in Refs. '13^14' reveal in fact, 
for the undriven case, that a perturbative approach in 
Jeff breaks down for strong qubit-detector coupling, and 
when the qubit and detector frequencies are comparable. 
Dephasing times at zero temperature have been deter- 
mined for the undriven spin-boson model with a struc- 
tured environment in Ref. Ila within a numerical flow 
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equation method. 

The interplay between the external driving and the 
dynamics of the coupled TSS-HO system yields to addi- 
tional multi-photon transitions, which can be explained 
only by considering the spectrum of the coupled sys- 
tem. These resonances have recently been observed 
experimentally^. If the time-scale of the HO does not 
play a role, the multi-photon resonances occur in the 
driven qubit solely, which happens when the driving fre- 
quency (or integer multiples of it) matches the character- 
istic energy scales of the qubit^ . Such multiphoton res- 
onances can be experimentally detected in an ac-driven 
flux qubit by measuring the asymptotic occupation prob- 
abilities of the qubit, as the dc-field is varied^^. In these 
experiments the resonances were obtained by matching 
the frequency of the ac-field with the qubit energy levels 
only, where the detector energy levels did not play a role. 
These qubit resonances, which have also been theoreti- 
cally investigated within a Bloch equation formalism in 
Ref. IITL could be explained in terms of intrinsic transi- 
tions in a driven spin-boson system with an unstructured 
environment. 

In this paper, we provide a comprehensive theoreti- 
cal description of the driven spin-boson system in the 
presence of a structured environment with one localized 
mode. Upon making use of the equivalence of this generic 
model with the model of a driven TSS coupled to an 
Ohmically damped HO, we first consider the experimen- 
tally most interesting case of low temperature and weak 
damping of the HO while the couphng between the TSS 
and the HO is kept arbitrary. In this regime, a Floquet- 
Born-Markov master equation can be established for the 
driven TSS-HO system. A restriction to the most rele- 
vant energy states allows the analytic calculation of the 
asymptotic TSS time-averaged population Poo, including 
the explicit shape of the resonance peaks and dips. We 
furthermore consider the case of strong damping and/or 
high temperature which is the complementary parame- 
ter regime. An analytic real-time path-integral approach 
within the non-interacting blip approximation for the 
driven TSS with the Lorentzian-shaped spectral density 
allows to analytically determine Poo as well. We compare 
the results obtained from closed analytic expressions with 
those of numerically exact real-time QUAPI calculations 
in both parameter regimes and find a very good agree- 
ment validating our analytical approaches. Finally, we 
consider the weakly damped TSS with the localized mode 
in the limit of large HO frequencies. Then, the localized 
mode acts as a high-frequency cutoff and the usual Ohmi- 
cally damped driven TSS is recovered. For this case, we 
employ an approximation valid for large driving frequen- 
cies and obtain a simple expression for the resonance line 
shapes for multi-photon transitions. Most importantly, 
we find that the width of the n-photon resonance scales 
with the n-th ordinary Bessel function. Parts of our re- 
sults have been published in a short work in Ref. 

The paper is organized as follows: In Sec.^ we present 
the theoretical model. Then, we treat the regime of weak 




FIG. 1: Schematic picture of the models we use. In (a) the 
TSS is coupled to an environment which has a peaked spectral 
density Joff(a;). In (6) the system is shown as a two-level 
system coupled to a harmonic oscillator which is itself coupled 
to an Ohmic environment with spectral density Johm('^) . 



damping and low temperatures in Sec. IIIII The comple- 
mentary regime of strong damping is investigated in Sec. 
HVl The subsequent Sec. contains the limit when the 
localized mode provides a high-frequency cut-off for the 
bath, and Sec. I VII the discussion of the results and the 
conclusions. Details of the specific evaluation of rate co- 
efficients are presented in Appendix ^ In Appendix ^ 
an expansion used in the strong coupling regime is elab- 
orated in detail. 



II. THE DRIVEN QUBIT COUPLED TO A 
MACROSCOPIC DETECTOR 



The driven TSS is described by the Hamiltonian 



Hgit) = — —ax 2~^'^^ ' 



(1) 



where tXi are Pauli matrices, HA is the tunnel splitting, 
and e{t) = Eq + scos(f2i) describes the combined effects 
of a time-dependent driving and the static bias Eq- bi 
the absence of ac-driving (s — 0), the level splitting of 
the isolated TSS is given by 



= hJel + A2 . 



(2) 



The detector can be associated as part of the TSS envi- 
ronment as a localized mode. This gives the spin-boson 
Hamiltonian Hsb reading2i2iiLiaiiiiiJi 

HsB{t) = HQ{t) + ]^a,h'^\k{bl + bk) + ^nwkblbk. 

k k 

(3) 



3 



Here bk and 6^ are annihilation and creation operators 
of the fc— th bath mode with frequency w^. The presence 
of the detector determines the shape of the spectral den- 
sity. Following Ref. 10, the dc-SQUID can be modeled 
as an effective inductance which is shunted with an on- 
chip capacitance. This gives rise to the effective spectral 
density 
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of the bath having a Lorentzian peak of width 7 = 
27rKf2p at the characteristic detector frequency flp. It 
behaves Ohmically at low frequencies with the dimen- 
sionless coupling strength a = limi^^o •^ff('^)/2w- The 
qubit dynamics is described by the reduced density op- 
erator p{t) obtained by tracing out the bath degrees 
of freedom. The relevant observable which corresponds 
to the experimentally measured switching probability 
of the SQUID bias current is the population difference 
P{t) := {az){t) = tr[p{t)az\ between the two localized 
states of the qubit. We focus on the asymptotic value 
averaged over one period of the external driving field, 
i.e.. Poo = \iTat-,oc{P{t))n- 

In the following, it will become clear that it is conve- 
nient to exploit the exact one-to-one mappingi^ of the 
Hamiltonian Q onto that of a driven TSS coupled to a 
single harmonic oscillator mode with frequency fJp with 
interaction strength g. The HO itself interacts with a set 
of mutually non-interacting harmonic oscillators. The 
corresponding total Hamiltonian is then 



HgoBit) — Hgoit) + Hqb 



(5) 



follows as g = ^py/a/JSK). Fig. ^ illustrates a sketch 
of the two equivalent descriptions of the system. Fig. 
^ shows the viewpoint where the localized mode is part 
of the environmental modes, while Fig. ^ depicts the 
perspective of the localized mode being part of the "cen- 
tral" quantum system which itself is coupled to an Ohmic 
environment. The equivalence of both standpoints has 
first been pointed out by Garg et ali^ in the context 
of electron transfer in chemical physics. As shown be- 
low, the first way is more convenient for the description 
in terms of analytic real-time path-integrals (Sec. IIV|1 . 
while the second viewpoint is more appropriate for the 
regime of weak-coupling and for the numeric treatment 
with QUAPI (see below). Note that the TSS reduced 
density operator p(t) is obtained after tracing out the 
degrees of freedom of the bath and of the HO. Further 
progress relies on approximations which depend on the 
choice of the various parameters. 



III. WEAK COUPLING: 
FLOQUET-BORN-MARKOV MASTER 
EQUATION 

If the coupling between the HO and the bath is weak, 
i.e., K ^ 1 , we can choose an approach in terms of a 
Born-Markov master equation in an extended Floquet 
description2i2SiSi. For a self-contained discussion, we 
shortly introduce below the required formalism of the 
Floquet theory. The interested reader can find more de- 
tails in the review in Ref. 0- 



A. Floquet formalism and 
Floquet-Born-Markovian master equation 
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HQo{t) 
Hob 



Hgit) + hga,{B^ + B) 
(St +B)Y, fi'^kibl + bk) + hwuhlbk 
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(6) 



where we have omitted the zero-point constant energy 
terms. Here, B and B^ are the annihilation and creation 
operators of the localized HO mode, while bk and are 
the corresponding bath mode operators. The spectral 
density of the continuous bath modes is now Ohmic with 
dimensionless damping strength n, i.e.. 



(7) 



where we have introduced a high-frequency Drude cut- 
off at frequency lvd. If is larger than all other energy 
scales the particular choice of cut-off does not influence 
the results at long times. The relation between g and a 



For systems with periodic driving it is convenient to 
use the Floquet formalism that allows to treat periodic 
forces of arbitrary strength and frequency^". It is based 
on the fact that the eigenstates of a periodic Hamiltonian 
HQo{t) — Hqoit + 2ti/VI) are of the form 



\Mt)) 



= \4>a.{t + 2TT/n)) 



(8) 



with the Floquet states \4)a{t)) being periodic in time 
(as is the Hamiltonian) and are called the Floquet or 
quasi-energies. They can be obtained from the eigenvalue 
equation 



d \ 

Hqoit) ~ ih- \ |0„(t)) = e„|0a(t)). (9) 



If the quasi-energy Sa is an eigenvalue with Flo- 
quet state |0Q.(t)), so is -I- nh^ with Floquet state 
exp {inQ.t)\(j)a{t)) ■ Both Floquet states correspond to the 
same physical state. Because of their periodicity both the 
Floquet states and the Hamiltonian can be written as a 
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Fourier series, i.e., 

71 

HQoit) = ^i/5^exp(inm). (10) 

n 

Substituting these Fourier decompositions in the eigen- 
value equation ^ givesSS 



(11) 



This aUows us to define the Floquet Hamiltonian Tiqa = 
Hqoit) — ifi''§i in niatrix form with the matrix elements 

{an\'HQo\bm) = (^fg'o^'^^jah + nhQSabSnm- (12) 

In the notation \an) , a refers to a basis in which to express 
the Hamiltonian Hgoit), while n refers to the Fourier 
coefficient. The eigenvectors of Hqo are the coefficients 

The dynamics of the system coupled to a harmonic 
bath is conveniently described by an equation of mo- 
tion for the density matrix p. Driving effects can be 
captured in an elegant way by formulating the equation 
of motion in the basis of Floquet states defined in Eq. 
(|5]) . For weak coupling to the environment, it is suf- 
ficient to include dissipative effects to lowest order in 
K. Within this approximation, a Floquet-Born-Markov 
master equation has been established^i^S'Si. We aver- 
age the 27r/r2-periodic coefficients of the master equation 
over one period of the driving, assuming that dissipative 
effects are relevant only on timescales much larger than 
27r/r2. This yields equations of motion for the reduced 
density matrix Papit) = {4'a{t)\p{t)\4'f3{t)) of the form 



Papit) = ^■^('^a ~ £l3)Pafsit) + X! ^aP.a' 0' Pa' 0' (t) 

a' 13' 



with the dissipative transition rates 

Laf3,a'l3' = ^^^(^aa'.ri + ^;3/3' ,Ti )^aa' ,n^/?'/3,-n 
n 

~^Pf3' ^ Np»a' ^nXaf3" ,-nXi3"a' ,n 
P" .n 

— Saa' ^ Na"i3' ,nXp'a" ,-nXa" f3,n- (14) 
a" ,71 

Here, we have defined 

XaP.n = 5:(0W|i? + St|0f+")), 



(15) 
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FIG. 2; Left: Quasi-energy spectrum ec,,k of the driven 
TSS-I-HO system vs dc-bias eo (in units of A). The quasi- 
energies are defined up to an integer multiple of hQ, i.e., 
£a,k ~ £a + khfl. Inset: Zoom of an anti-crossing. Right: 
Poo exhibits resonance dips corresponding to quasi-energy 
level anti-crossings. Parameters are Q = lOA, s = 4:A,g = 
0.4A, Q.p = 4A, K = 0.014 and fcsT = O.lfiA. 



We have neglected the weak quasi-energy shifts, which 
are of first order in the coupling to the environment. In 
the sequel, we will see from a comparison with exact nu- 
merical results that this approximation is well justified. 
In order to be able to solve Eq. l(T^ . it is necessary to de- 
termine the Floquet quasi-energies and Floquet states 
How they can be determined perturbatively, is 
shown in the following subsection. 

B. Van Vleck perturbation theory 

First, we have to specify the basis for the Floquet 
Hamiltonian according to Eq. For the TSS-I-HO 

Hamiltonian Hqq, we use the basis \an) defined via the 
single particle product state \a) — \g/em) with \g/e) be- 
ing the ground/excited state of the qubit, \m) the HO 
eigenstate, and n the corresponding Fourier index. In 
detail, this implies that we can divide the Hamiltonian 
into a diagonal part 

0-iQo)gmn,gmn = h[-iy /2 + Uinp + nfl] , 
{'HQo)emn,emn = h[u / 2 + mflp + nfl] , (16) 

and an off-diagonal part 

{'HQo)an,bk = Van,bk, for fl 7^ 6, 71 7^ /c , (17) 

which has non-zero elements. They read 



g{e)ln,g{e)mn 



+(-) 



{\/m+ 1(5;, m+i + \/r+T(5;+i,m)%£o 



g(e)ln.e(g)mn 



(Vm + l(5;,TO+i + + l(S;+i,m)%A 



g(e)7nn.g(e)mk 



-(+■ 



g(e)mn.e{g)mk 



(18) 



In the remainder of this section we will assume that the 
elements of V are small compared to the diagonal el- 
ements of Ti-QOi which is justified if the coupling g be- 
tween TSS and HO and the driving amplitude s are small 
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compared to the other energy scales, i.e., s,g <^ D,, v, Qp. 
This is the case in reahstic experimental device a^'^^ . The 
Fourier index n ranges from — cx) to cx) and m from to 
oo. The eigenvalues of the Floquet Hamiltonian following 
from Eqs. Ijl?!) and l|18() have to be calculated numerically 
for a particular cut-off rimax and rrimax- In Fig. |21 the 
numerically obtained quasi-energy spectrum is shown as 
a function of the static bias £q for the case Wmax = 4 and 
I'^maxl = 8. We find that for some values of the bias eg 
avoided crossings of the quasi-energy levels occur when 
two diagonal elements of Ti-qo have approximately the 
same values, i.e., when the condition 



Ea 



{T^Qo)an,an — (Ji-Qo)bm,bm = + 0{V^) (19) 



is fulfilled. It follows from Eq. H17() that this happens 
when at least one of the two conditions 



(20) 



is fulfilled. At these avoided crossings the Floquet spec- 
trum has quasi-degeneracies and as a consequence there 
are transitions between the different Floquet states. As 
it turns out below, this results in resonant peaks/dips 
in the stationary averaged population difference Poo , cf. 
Fig.H 

Since we are interested in describing the resonance line 
shape for Poo, we have to determine the quasi-energies 
and Floquet states around a resonance, i.e., around an 
avoided crossing. For this, we use an approach which 
is perturbative in V . The unperturbed Hamiltonian is 
diagonal and, close to an avoided crossing, (nearly) de- 
generate. An appropriate perturbative method is the 
Van Vleck perturbation theorji^^ suitable for Hamilto- 
nians for which the unperturbed spectrum has groups 
of (nearly) degenerate eigenvalues, well separated in en- 
ergy space. An example of such a spectrum is shown 
in Fig. 121 This method defines a unitary transforma- 
tion which transforms the Hamiltonian into an effective 
block-diagonal one. The effective Hamiltonian then has 
the same eigenvalues as the original Hamiltonian, with 
the quasi-degenerate eigenvalues in one common block. 

The effective Hamiltonian can be written as 



Tic 



e'^HQoe 



■iS 



(21) 



In Ref. it is shown how to obtain S systematically for 
every order in the perturbation. The small parameter is 
y/Ai?2 (see Fig. O. Eigenvalues within one block can 
be arbitrarily close. This means that we can also use the 
expansion at resonance. We derive the expressions up 
to the second order in the perturbation. Two different 
cases are relevant: For the case when {'HQo)an.an and 



AE, 



FIG. 3: Typical energy spectrum suited for the Van Vleck 
perturbation theory: Different groups of (nearly) degener- 
ate levels of eigenenergies are well separated in energy (i.e., 
AEi < AEi). 



iJiQo)bm,bm are not nearly degenerate we find 



. (-f(l) * an^bra 



a7i,bm 
(2) 



Ean,bm 

Van,ckVck,bn 

c.k 



2Ebn 



Ec 



Ec 



k.bn 



(22) 



where the superscript indicates the order of perturba- 
tion theory. For the second case when {'HQo)an,an and 



(1) 



{'HQo)bm,bm are nearly degenerate, we find iS^ 

^'-'an,bm — 

In turn, the matrix elements of the n-th order term 

(n) 

n^ff of the effective Hamiltonian can be calculated, 
again for both cases. When (TiQo)an.,an and {TiQo)bm,bm 
arc not nearly degenerate, we find ('Woff)|,^\m = 

(^cfr)L fcm = 0. For the second case, when (Hqo) an, an 
and {'HQo)bm.bm are nearly degenerate, one finds 



(Wcff) 



(1) 

an^b7n 



(^off)f2 



bm 



an.brn i 



E, 



an.ck 



E, 



(23) 



Since the effective block-diagonal Hamiltonian conse- 
quently only has 2x2 blocks, it is easy to diagonalize 
it. To obtain the eigenvectors of the original Floquet 
Hamiltonian (see Eqs. H17ll8(l V the inverse of the trans- 
formation defined in Eq. H21|l has to be performed on the 
eigenvectors. There is an infinite number of quasi-energy 
levels and Floquet states. However, because the eigenval- 
ues £q, and Ea + nhfl represent the same physical state, 
only one of them has to be considered. Still there is an 
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infinite number of levels because the Hilbert space of the 
HO Hamiltonian is infinite dimensional. Nevertheless, for 
practical calculations, only the relevant HO levels have 
to be taken into account. When there is a resonance 
between the states |e/g,0,rt) and je/g, Z,n+ fc), then at 
least the first I levels of the HO play a role. Higher levels 
can be omitted if one is interested in low temperatures, 
which is commonly the case, since for low temperatures 
their occupation number will be very small. 

C. Line shape of the resonant peak/dip 

To obtain the line shape of the resonant peak/dip in 
Poo , we have to determine the stationary solution of Eq. 
(|13|l . Depending on the number Umax of Floquet states 
taken into consideration, this might be considerably dif- 
ficult. One facilitation might arise due to symmetries, 
i.e., elements oi Lap^a'fj' being of the same size. Another 
possibility appropriate at low temperatures might be to 
neglect some of the dissipative transition rates. More- 
over, a further possible approximation can be applied for 
Pa/sioo), if Ea and e/3 are not nearly-degenerate eigen- 
values. In that case Eq. — in Eq. H13|l is much larger 
than the coefficients Lapa'p', since the coupling to the 
Ohmic environment is assumed to be weak. This, in 
turn, allows to make the partial secular approximation 
by setting papioo) = 0. After the reduced density ma- 
trix in the Floquet basis is known, it is straightforward 



to calculate Pc 



D. Example: The first blue sideband 

As an example we will derive an analytical expression 
for the resonant dip a,t v ^ — Q,p which is called the 
first blue sideband. For this case, the matrix elements 
{'H.Qo)gOn+i,gOn+i and {HQo)ein,ein are nearly degener- 
ate, i.e., 

-1^/2+ {n+l)nK 1^/2 + Qp + nQ. (24) 

As the resonance occurs between two states which differ 
only by one oscillator quantum, we only take into account 
one excited level of the oscillator. We expect that this 
is a reasonable approximation for not too strong driving 
and low temperatures. The validity of this approximation 
will be checked against numerically exact results in the 
end. The elements of the transformation matrix S follow 
as SgOn+i,ein = Sein,gOn+i = 0, whilc thc remaining 
elements can be calculated straightforwardly using Eq. 
(|22|l and they will not be given here explicitly. 

Since we include one HO excited energy 
level, we have four physically different eigen- 
states. Hence, we can express Ti.Qo in the basis 
{|e,l,-l),|5,0,0),|(?,l,0),|e,0,0)}. Performing the 
transformation defined in Eq. H21|) . we obtain the 
effective Hamiltonian in this basis as 



Ai 







Wi 



-Ai 



W2 











Ws 



(25) 



W4 



The matrix elements are calculated using Eq. H23|l . They 
read 



Aeogs [n^ + »^ + i^j-H + ftp)] 

4i^{n - iy)nnp{i' + Up) 



W3 = -Wi = 





h 




ly^ flp 


iy^{v + ap) 




^0.9 


A2.g2 


A2s2 


v'^ flp 


h''^(iy — ilp) 





(26) 



The eigenvalues of the Hamiltonian H25|) are the relevant 
quasi-energies, and they are readily obtained by diago- 



nalization as 

-j---,+W4. (27) 

From these formulas it follows that 6 = i'—rt+il,p+2Wi is 
a measure of how far the system is off resonance. For S — 
0, the quasi-energies ei and £2 show an avoided crossing 
of size hAi. Note that Eq. implies that Ai is the 
Rabi frequency at the blue sideband. 

The eigenvectors, which are the Floquet states, of the 
4x4 effective block-diagonal matrix in Eq. (|25|l are easily 
obtained by performing the corresponding inverse trans- 
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formation. We find, withtan^ = 2|Ai|/5, the eigenstates 

= e-*^[sin(^/2)e-^"*|e,l)+ cos (0/2)15,0)], 

102) = e-*^[cos(0/2)e-'^*|e,l)-sin((?/2)|5,O)], 

l'/'3/4) = e-*^|5/e,l/0). (28) 

We have used the inverse transformation of Eq. 1)10(1 to 
illustrate the time-dependence explicitly. Next, we cal- 
culate the rates given in Eq. H14() up to second order in 
V. 

The quasi-energies ei and £2 are quasi-degenerate. To 
be definite, we assume a partial secular approximation: 
We set almost all off-diagonal elements of p to zero but 
keep Pi2(oo) and /32i(oo) = p\2{oo) different from zero. 
This allows to simplify the master equation ((13(1 . The 
stationary solutions are determined by the conditions 



= 



/3 

i 



Lota, PPPppioo) + {Laa,12 + iaQ,21 )Rc [/3i2 (oo)] , 



= --^(£1 - £2)^12(00) + Li2^aaPaa{oo) 



+-^12,12^12(00) -I- L12, 21^^2(00). 



(29) 



It is most convenient to use the symmetry properties of 
the corresponding rates which are specified for this par- 
ticular example in the Eqs. ((Aip in Appendix 1X1 In turn, 
there are eight independent rates associated to all pos- 
sible transitions. They are explicitly given in the Eqs. 
(IX2I1 . 

First we consider the rates exactly at resonance (5 = 0. 
Since then sin^(0/2) = cos^(0/2) = 1/2, all rates contain 
a term which is of zeroth order in y. If we neglect the 
small second order terms, we find 

^ -^33, 21 



^44, 21 



^22, 44 


— -^11,44 — 


^33, 22 


— ^33,11 




= -^21,44 = 


N{hnp) , 


-^22,33 


= -^11,33 = 


-^44,22 


= ^44, 11 




= ^^21, 33 


= N{^ 




ill, 21 


= -^22,21 = 


^21, 22 


= ^21,11 



1 



[N{nnp) - N{-nnp)] , 



ii242 = -Ni-hrip) - N{hnp) . (30) 

Solving Eq. (jSHl together with (jSOJ finally yields 

N{-hnp)N{hnp) 



Pll(oo) = P22(oo) 

P33(oo) = 
P44(oo) = 



N{hnp)^ 



[N{-nnp) + N{hnp)]^ ' 



[Ni^hQp) + N{hnp)]^ ' 



pi2(oo) = 



(31) 



Eventually, this gives the simple result at resonance, d = 

0, 



which implies a complete inversion of population at low 
temperatures! We will discuss the physics of this in sec- 
tion ^^E| Note that no further assumption on the tem- 
perature was made while deriving this formula. 

Next we will derive an expression for the peak shape 
around the resonance. For this, we assume low tem- 
peratures, i.e., ksT/h <C Qp^fl^v. This allows us to 
set Nikflp) = N{hn) = N{hv) = 0. Far enough 
away from resonance, it is appropriate to assume that 
/9i2(oo) « /92i(oo) « 0, and sin(6'/2) « 9/2. Thus, it fol- 
lows from Eq. ((Aip that there are only four independent 
rates in this case, namely £44,22, ^22, 44, ^44, 11 and £11^44. 
Within our approximations, we find that ^22,44 = 
So only three rates are relevant which read 



^44,22 = 2cos2(6l/2) « 2, 
Lii,44 = 2L^{eiAfi)cos^{9/2) 



8A-'g'n-'p 



^qV^i,4,uy ^^'J , o o2\2 ' 

[V^ — V\lp)'^ 

L44,ii = 2sin2(6'/2) w 6(2/2 « 2A?/(5^ (33) 

where the quantity is given in Appendix ^ Note 
that L44 22 3> ^44^11, L11.44. In this limit we find for the 
asymptotic density matrix elements 



pii(oo) = 



Li 



1,44 



/022(oo) 
P44(oo) 

which gives the central result 
£0 -^11.44 — Lu ll 



ill, 44 + -^44,11 

P33(oo) = 0, 

1 - pii(oo) , 



(34) 



Poo = 




44,11 



2A2i.2(,,2 



(35) 



A comparison between the result of this formula and dif- 
ferent numerical results, including those of an exact nu- 
merical ab-initio real-time QUAPI calculationSiiS42&, is 
shown in Fig. 21 For the QUAPI-simulations, we have 
used the optimized parameters2& At = 0.23/ A, M = 12 
and K = 1. Moreover, we have applied an exponen- 
tial cut-off for the Ohmic bath with a cut-off frequency 
LUc — lOA (since we are considering long-time station- 
ary results, the explicit shape of the cut-off is irrelevant). 
Note that the picture of the TSS+HO being the central 
quantum system which is coupled to an Ohmic environ- 
ment is particularly suited for QUAPI since the coherent 
dynamics of the central quantum system is treated ex- 
actly. A very good agreement, even near resonance, is 
found among all the used numerical schemes. 



E. Results and discussion 

To get a qualitative understanding of the results it is 
important to realize how the rates Lap^a'P' are formed. 
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All Floquet states are superpositions of the four unper- 
turbed states \g,0), Iff, 1), |e,0), |e, 1). The rates Lose, 
Lq, iq,osc defined in Eq. 1A2|I describe the time scales of 
the transitions between the different unperturbed states. 
In this discussion these will be called the basic transition 
rates. If we want to calculate the rates Laa'^pp' which 
describe transitions between different Floquet states, we 
have to multiply the basic rates by the square of the am- 
plitudes in the superpositions. For example the Floquet 
state consists of an unperturbed state |e,l) with 
amplitude sin(6'/2), and |(?, 0) with amplitude cos(0/2). 
The rate ^44,11 describing the dissipative transition from 
|0i) to |(^4) has a term which is the basic transition rate 
from |e, 1) to |e,0) (Lose) multiphed by sin^(0/2), and 
a term which is the product of the basic transition rate 
from Iff, 0) to |e, 0) (Lq) and cos^(0/2). In this quahtative 
explanation we can neglect the amplitudes of the other 
states in |(/)i) and |04) which are 0{V'^). 

Now we consider the rates at resonance. Here, the am- 
plitudes sin(0/2) and cos(0/2) are equal and we only have 
to consider the largest basic rate which is Lose- The domi- 
nant transition from |0i) to |04) is via the basic transition 
from |e, 1) to |e,0) and it is fast (both amplitude in the 
superposition and rate Lose are of order one). The same 
holds for the transitions between |02) to \4>4). There will 
also be a transition between and \4'i), |02) via the 
basic transition between |ff, 1) and Iff, 0). Since only the 
decay of the oscillator plays a role, the stationary state 
is a thermally equilibrated mixture between |e/ff, 1) and 
|e/ff,0) which is described by Eq. 132|) . 

Away from resonance, the amplitude sin(0/2) becomes 
small {'^ Ai) implying that other basic transitions start 
to play a role. The transition from |(/)2) to |(/)4) is still 
dominated by the basic transition from |e, 1) to |e, 0) with 
large amplitude cos(0/2). Within our assumption of low 
temperatures used to derive Eq. (|^ . the state |02) is 
weakly populated after long times. The same holds for 
\4>3)- Between \<j)i) and |04) two basic transitions are im- 
portant. The first is a transition from |e, 1) to |e,0) with 
large rate (i.e., fast) but low amplitude (sin(6'/2)) (i.e., 
rare). For low temperatures, this rate describes transi- 
tions from |0i) to \(j)4) and it dominates L^^^n in Eq. 
(|33|l . The other basic transition occurs via the unper- 
turbed states |e, 0) and |ff, 0). The basic rate Lq is small, 
but it has a large amplitude (cos(^/2)). For low tem- 
peratures this mechanism causes transitions from |04) to 
|0i) and it is the dominant part of Lii^44. The stationary 
state is described by the ratio of the rates as described 
by Eq. (|35|l . Note that since both rates scale with 
the final result does not depend on g, for small g. For 
other resonances we find a different eigenvalue splitting 
and the peak shape will depend on g. 

We will now try to give a more physical insight in the 
nature of this resonance. First consider the case when 
T = and we are exactly on resonance. The driving 
induces transitions from |ff,0) to |e, 1) while the direct 
coupling of the HO to the environment will cause a fast 
decay of the population from |e, 1) to |e,0). This tran- 
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FIG. 4: Poo vs eo (in units of A) around the peak at = — 
r^p. The solid lines are the analytical prediction 13511 for (a) 
g = 0.05A, (6) g = 0.2A, (c) g = 0.4A. The triangles are the 
results of a Floquet-Bloch-Redfield simulation, cf. Eq. 
with one (upward triangles) and two (downward triangles) 
HO levels taken into account. The circles in (&) are the results 
from a QUAPI simulation with six HO levels (see text). We 
choose s = 2A, n = lOA, k = 0.014, fcsT = O.lftA. 



sition from |ff, 0) to |e, 0) via driving and decay has to 
compete with the decay from |e, 0) to |ff, 0), but the last 
process is much slower because the TSS is not directly 
coupled to the environment. So all the population is in 
|e, 0) and there is a complete inversion of population. 
For T 7^ there will be a thermal equilibrium between 
ground and excited state of the oscillator. 

When the system is not exactly at resonance the driv- 
ing induced transitions are much slower oc g^ and the 
decay of the oscillator is still fast. This means that the 
transition from |ff , 0) to |e, 0) is slower than at resonance. 
The time associated with the decay from |e, 0) to |ff,0) 
is also cx ff^ and the ratio of the time scales of the two 
processes gives the ratio of the populations of |e, 0) and 
|ff,0) (at T = 0, for higher T the states |ff/e, 1) are also 
populated) . This ratio is independent of g and so is Poo ■ 

A similar analysis can be performed for the first red 
sideband at = -I- f2p. At resonance it yields Poo = 
a tanh(^f^) + 0{V'^), which is very close to thermal 
equilibrium for low T. 

For the resonance at = Sip, only the oscillator is ex- 
cited. After having traced it out, we expect just thermal 
equilibrium given by Poo — ^ tanh( 2^^^ ) ■ 



IV. STRONG COUPLING: NIBA 

In the complementary regime of large environmental 
coupling and/or high temperatures it is convenient to 
employ model (a) , and it is appropriate to treat the sys- 
tem's dynamics within the nonintcracting-blip approxi- 
mation (NIBA)^. The NIBA is non-perturbative in the 
coupling a but perturbative in the tunneling splitting A. 
It is a good approximation for sufficiently high tempera- 
tures and/or dissipative strength, or for symmetric sys- 
tems. Within the NIBA and in the limit of large driving 
frequencies $7 A, one finds^ 



Poc 



MO) 

fco+(0) 



(36) 
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Here, 

V(0) 
k+iO) 



NT /'2s . 
dtn (i) sm (£oi)>^o I — sm 







/ dth+{t) cos {eot)jQ 
'o 



n 



(37) 



with Jo being the zcroth order Bessel function. Dissipa- 
tive effects of the environment are captured by the terms 



h+{t) = e-«'(*)cos[Q"(t)] 
h-{t) = e-'J'W sin[0"(i)] . 



(38) 



Here, Q'{t) and Q"{t) are the real and imaginary parts 
of the bath correlation function 



Q{t) = 



J{ijj) cosh(w/3/2) - cosh[w(/3/2 - it)] 



70 w2 sinh(w/3/2) 

For the peaked spectral density given in Eq. one finds 

Q\t) = Q;(t)-e-"[riC0s(f2pt)+r2sin(npt)], 
Q"{t) ^ Ai- e-^^[Ai cosiClpt) + A2 smiClpt)] . (39) 

Here, f3 = h/ksT, T = ttk^p, 0.1 = ^^- and 

'smh{pnp)t sin(/3r)i 



Q'lit) = 



4a 



2 00 



2cr 



1] +i 



(40) 



where i/„ = 2TTn/l3. Moreover, C = cosh(/3f2p) — cos(/3r), 
CYi/2 = T^2/isinh(/3f]p)-Ai/2sin(/3r), = a^(r2- 
f]^)/(2rf]p), ^1 = 7ra. As follows from Eq. Q' and 
Q" display damped oscillations with frequency Clp (cf. 
Fig. \^) which are not present for a pure Ohmic spec- 
trum. It is the interplay between these oscillations and 
the driving field which induces the extra resonances in 
p 

To proceed, we rewrite the kernels (0) in a more con- 
venient form. In the integrand of Eq. H37|l the functions 
cos[Q"(<)], sin[g"(i)] and e-<3'(*)+'3i(*) oscillate with fre- 
quency flp and we can expand them as 

oc 

cos[Q"{t)]~ ^ [D„iCOs{mapt) + E„iSm{mrtpt)] , 

m— — 00 
00 

sin[(5"(t)] = ^ [F„iCos{mdpt) + GmSin{mClpt)] , 

m— — 00 
00 

g~Q'(t)+Q'i(«) = [H^cos{mnpt) + K,nSin{mnpt)] . 



0.5 
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FIG. 5: Poo vs eo (in units of A). The solid line is the NIBA 
prediction, while the circles are from a QUAPI simulation 
with 6 HO levels (5 = 3A, s = 4A, Q = lOA, k. = 0.014, 
ksT = 0.5fiA, rip = 4A). Inset (a): NIBA result for ksT = 
2hA. The arrows indicates the first red sideband at u = 
^ + Qp. Inset (b): Q'{t) vs t shows damped oscillations. 



these expansions into Eq. H37|l . and also using the Fourier 
representation of Jq (|f sin^), we find 



oc 00 



C30 



fco"(0)= 1^ 1^ / rfte-'JiW/ijt). (42) 

m— — 00 n— — 00 







Here, Smn = £o — mQp — nQ, and 

Imnit) = Ktm(^)cos(e™„i)±4,„(i)sin(e,„„t)] 
= Jn (^) J™(e-"c.i)cos(m</>)(-*)™e-^^\ 



c- = 72 



(£) J™(e-"c.i)sin(m(/))(-z)"e-^^^ , (43) 
with Jn being a Bessel function of order n, and 



1 = ^{A^-iY^Y + {A2^iY2Y. 
A2 - iY2 



tan (h = — 



Ai - iYi 



(44) 



Thus, from Eq. H43|l . we expect resonances when e„,„ — 0. 

In the limit r/57p ^ 1 and for not too large T (i.e., 
cos(/3r) ^ cosh(/3ilp)), we find that 



tan(m(5 



i tanh 



m(3Qp 



Inserting this into Eq. (|43|l . we obtain 

.tanhril^^c;t.„=c„„. 



(41) 



The coefficients Dm, Em, Fm, Gm, Hm and Km are time 
dependent and they are given in Appendix ^ Inserting 



(45) 



(46) 



If the environmental mode is enough localized (i.e., the 
integrand of Eq. (|37|l is only damped after several oscil- 
lations), we expect that the sum in Eq. (|42|l is dominated 
by the coefficient of cos Enmt if £mn = 0. This means that 

/+„(i) « Re[c+„(i)] 

fmnit) « tanh (ai^) Rc[c+„(t)] , (47) 
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which leads to 



tanh 



mj3VLp 



(48) 



Without driving we only have terms with n = and 
£mO = implies that eg = mfip. In that case Eq. 
gives the NIB A thermal equilibrium value. Hence, in or- 
der to find resonances we need to apply driving. For 
"conventional" resonances at — '"^^i we put m to 
zero and we find Poo ~ 0, as predicted for unstructured 
environmentsiii^. Finally, for eg = nO. ± mllp, we re- 
cover Poo ~ ± tanh(TO/?rip/2), as was also found within 
the Floquet-Born-Markov approach, cf. (|32|l . Results of a 
numerical evaluation of Pqo are shown in Fig.jS] using the 
NIB A result (|43f) , as well as the exact ab-initio real-time 
QUAPI methodSiiS^. Resonance dips are observed at 
£q = £f) = — VLp and Eq = ^ ~ 2D,p. For kgT ~ Mlp^ 



we also find the first red sideband dX £q = VL 
inset (a). 



see 



In the limit when the frequency ilp of the HO is much 
larger than the effective TSS level splitting the peak in 
the spectral density at Op acts as a high-frequency cut-off 
for an effective Ohmic bath, see Eq. I^J. In other words, 
the oscillations in the correlation functions, see Sec. IIVI 
which occur on a time-scale are very fast and can 
be averaged out when only the long-time dynamics is 
of interest and short-time effects are not considered. In 
this limit, the standard driven and Ohmically damped 
spin-boson model2*2i is recovered. In the regime of weak 
damping a <C 1, the stationary population difference Poo 



has been determined within the assumption of large driv- 
ing frequencies f2 ^ A, ... upon using analytic real-time 
path integral methods in Ref . lU In this Section, we use 
this high-frequency approximation of Ref. '2^ as a starting 
point, and derive a closed simple analytic expression for 
the peak shape of the "common" multi-photon resonance. 
Most importantly, we find the scaling of the width of the 
n-photon resonance as the n-th Bessel function J„(s/r2). 
This scaling behavior has been observed experimentally 
in superconducting flux qubit devicesi^. 

The central issue in finding a closed analytic expression 
for Poo is to find the roots z9„ of the pole equationSi 



+ 00 



-t-oo 



+ 00 



m— — oo,7ri^n 



(49) 

Here, e„ — £q — nfl is the photon-induced bias and A„ = 
I J„(s/0)|A is the field-dressed tunneling splitting of the 
TSS, where Jn{x) is the n-th ordinary Bessel function. 
Considering the n-photon resonance, we numerically find 
that, up to extremely high numerical precision, the roots 
of Eq. H49|l are given by 



^0 



^0 ' 



Je^ + A2 



(50) 



Plugging Eqs. IjSOl) in the expressions for Poo given in 
Ref. I27L see Eqs. (6) and (7) therein, we find a closed 
expression for the lineshape of the n-photon resonance 
to be 



,2 \ ■/Ohm('?0) 



1?0 



|A2(e2 

I nV^n 



-0) 



1?„ 



.2 2 



;)'^Ohm(^^o)coth 



2kBT 



A2e2(,2+A2_,2)j^^^^(^^ 



icoth^ 



(51) 



This result can be simplified upon observing that the sec- 
ond term in the numerator is small if the driving is not 
too large, since then A^ <C Ag. Moreover, we are in- 
terested in the regime eg 3> A which is the saturation 
regime implying that « eg and at low temperatures. 
Then, away from the resonance point at nO w Eq, the 
second term in the denominator in Eq. H51|l can be ne- 
glected and one recovers the standard result, if one uses 
that Aq ~ A which is fulfilled for weak driving. It reads 



: tanh 



(52) 



which gives the correct result away from any n-photon 
resonance. For the case at the n-photon resonance at 



nVL £0, one finds a Lorentzian line shape, i.e. 

A2(£o-nr!)2 



A2(eo - nnf + 2el/^lkBT/h 



(53) 



where we have expanded the second coth term in the 
denominator (see Eq. up to lowest order in the ar- 

gument, which is appropriate since "dn is small at reso- 
nance (and at low temperature). The linewidth of the 
Lorentzian peak can be calculated as the full width at 
half maximum (FWHM) 



Ae(") = 2 




An 

Ao 



(54) 
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Note that this result obtained from the high-frequency 
approximation is independent of the damping constant. 
Moreover, the leading term is the first term under the 
square root in Eq. I|54|) . Note furthermore that for the 
case of infinitesimal driving, the FWHM is not correctly 
reproduced by Eq. (|54|l since it would approach zero. 
However, as it is known from NMR within a treatment 
in terms of the Bloch equation, in this case, the FWHM 
is dominated by the dephasingii, i.e., 

A4Voch = V^^ + f^flVri^, (55) 

where — 7racoth(fi,z^/2fcBT)A^/i/ is the relaxation 
rate and F^ = TR/2 + 2TTa{el/i^^)kBT/tf. Both rates are 
of first order in the damping strength a. Moreover, fi/j 
is the (single-photon) Rabi frequency. Hence, we have to 
include the dephasing rate F^ in Eq. (|53|) , since it cannot 
be reproduced by our weak-coupling approach which is 
only of first order in a. This finally yields in leading 
order in the driving strength 

A.<">^2^rj+(f^)%„0^f ^ (56, 

As follows from Eq. ifS^ . the FWHM of the n-photon 
resonance scales with the n-th ordinary Bessel function, 
i.e., Ae^") ~ J„(s/f2) as also confirmed by experimental 
measurements^^. 



VI. CONCLUSIONS 

In conclusion we have investigated the problem of a 
quantum mechanical driven two-state system being cou- 
pled to a structured environment which has a localized 
mode at a frequency f2p but behaves Ohmically at low 
frequencies. We have studied two complementary pa- 
rameter regimes of weak and strong coupling to the envi- 
ronment. The interplay of the driving and the localized 
mode gives additional features like resonant peaks/dips 
in the asymptotic averaged population difference Poo • We 
have calculated analytically the lineshape of the reso- 
nances in various parameter regimes and have obtained 



simple closed expressions for the particular example of 
the first blue sideband. We also include the discussion of 
how the results are generalized for any sideband. More- 
over, we have elaborated the limit when the localized 
mode acts as a high-frequency cut-off. Then, the full 
width at half maximum of the n-photon resonance has 
been shown to scale with the n-th ordinary Bessel func- 
tion. 

Our model finds as well applications in the field of cav- 
ity quantum electrodynamics (CQED) with solid state 
structures'^. Most interestingly, the strong coupling 
limit of CQED could be reached in superconducting elec- 
trical circuits, with perspective applications ahead. 

Finally, we note that a related experiment has been 
reported recently by Wallraff et alM>. There, a qubit was 
realized in the form of a Cooper pair box which couples 
to a single mode of a cavity which is damped. The prop- 
erties of the TSS-HO were probed spectroscopically by 
measuring the transmission of the resonator. In other 
words, a driven HO was considered while the TSS was 
kept static. In contrast to that system, here, the TSS 
was time-dependent while the HO is treated as static. 
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APPENDIX A: SYMMETRY PROPERTIES FOR 
THE DISSIPATIVE RATES FOR THE FIRST 
BLUE SIDEBAND 

In order to evaluate the stationary averaged popula- 
tion difference Poo, the rate coefficients La/3,a'p' have to 
be determined explicitly. For the example of the reso- 
nance at ~ fl — flp (first blue sideband) considered 
in this work, we find that the rate coefficients fulfill the 
symmetry relations 



Pll,22 


— P224I — 


i33,44 — 


i44,33 = 0, ill, 11 = — (i33,ll + i44,ll), 


p22,22 


= ^(-^33,22 


+ i44,22 


), i33,33 = — (ill,33 + i44,33). 


^44, 44 


= ^(ill,44 


+ i22,44 


), il2,12 = (i22,22 + ill,ll)/2, ill, 33 = i44,22 


^22, 33 


= ^44, 11, 


i33,22 — 


ill, 44, i33,ll — i22,44, i22,21 — i21,ll. 


ill, 21 


= ^21, 22, 


i44,21 = 


— i21,33, i33,21 = — i21,44, i22,12 = i22,21. 


ill, 12 


= ill, 21, 


i44,12 = 


i44,21, i33,21 = i33,12, il2,21 = 0. 



r 



As a consequence, there are eight independent rates given 
by 
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-^44,22 — 2Losc(e4,2,-l) COS^ 



2-Z^q(e4,2,o) sin^ 



^22,44 = 2Losc (£2,4,1) COS ( - ) + 2Lq(e2,4,o) sin"- 



-^44,11 
-^11,44 

-^21,22 
^21,11 
-^21,44 

^21, 33 

with 



2Losc(e4,i,-i)sin^ ( ^ ) + 2Lq(£4,i,o) cos^ 



- Lq,osc(£4,2,-l)sin6', 

iq,osc (£2,4,1) sin6' , 
'T' -i^q osc (£4,1,-1) sin 6' , 
2iosc(£i,4,i) sin^ l^^j + 2iq(ei,4^o)cos^ j + -Lq^osc(£i,4,i) sin6' , 

1 (Losc(£3,2,o) - iq(£3,2,-l) - -^osc (£4,2,-1) + iq(£4,2,o)) sill 6* + ^ (Lq,osc (£3,2,o) - iq,osc(£4,2,-l)) COS 6* , 
1 

2 
1 

2 

2 (■^q(£l,3,l) +iq(£2,3,l) - ^osc(£l,3,0) " ^osc (£2,3,o) ) sill 6* - - (iq,osc (£l,3,o) + iq,osc (£2,3,o) ) COS ( 



^ (Losc(£3,l,0) - iq(£3,l,-l) - -Lose (£4,1,-1) + iq(£4,l,o)) sin 6* + ^(iq,osc(£3,l,0) " iq,osc(£4,l,-l)) COS 6* , 
^ (Losc(£l,4,l) + iosc(£2,4,l) - -Lq(£l,4,o) " -Lq(£2,4,o)) sill 6* + ^ (iq,osc (£l,4, 1 ) + iq,osc (£2,4, 1 ) ) COS 6* , 



iq(£fcim) = (e, 1, 0|e'^Xe"** \g, 1, 0)^iV(£fc,„0 = (e, 0, Ole^^Xe"*^ \g, 0, OyN{ekim) = 



I/2(i/2_^)2 



Losciekim) = {g,l,0\e'^Xe-'^\g,0,0)^Niekir„) - {e,l,0\e'^ Xe-'^\e,0,OfN{skim) 

^ + 1 I A^(£fei™) , 



^ / 4g^(A^(i.^-20g)-(^^ -17g)^^ 

iq,osc(£fei™) = -2(e, 1, -l|e'^Xe-^^|g, 1, 0)(5, 1, 0\e'^ Xe-'^\g, 0, 0)iV(£fe,„) 
= 2(5, 0, l|e'^Xe-'^|e, 0, 0) (e, 0, 0\e'^ Xe-'^'\e, 1, 0)A^(£fez™) 

A£o.g5((» + + 2»g + + »p))jV(£fch„) 

~ 2i^(iy-fi)17(i/-f7-f^p)f^p(i^ + f7p) 

N{ekim) = iV(£fc - £; + mf7) . 



(A2) 



r 



For the expansion of cos[Q"{t)], we find 



Note that iosc is the rate containing the zeroth order 
term in g and s. It is related to the transition between 
two states differing by one oscillator quantum. This de- 
cay is of zeroth order (hence fast) because the oscillator is 
coupled directly to the environment. Moreover, Lq gives 
the rate for transitions between the excited and ground 
state of the TSS with the HO remaining in the same 
state, and iq,osc is related to the transition where both 

the qubit and'the oscillator exchange energy with the en- where we have introduced 
vironment. Note that this transition is induced by the 
driving and involves one photon. 



-D2m+i = (-l)'"sin(Ai)J2™+i(A)cos[(2m + l)X], 

D2m = (-l)'"cos(Ai)J2,„(A)cos(2mX), 

E2,n+i = (-l)'"sin(Ai)J2™+i(^)sin[(2TO+l)X], 

E2m = (-l)'"cos(Ai)J2„,(A)sin(2mX), (Bl) 



A 



e-''JAl+Al, 



APPENDIX B: COEFFICIENTS FOR THE 
KERNELS k^{0) 



sinX = A^/^Ai + Ai. 
In the same way, the expansion of sin[Q"(i)] gives 



(B2) 



In Sec. IIVI we have introduced an expansion of the os- 
cillating functions given in Eq. imil. In this appendix we 
summarize the corresponding coefficients for complete- 
ness. 



F2„+i = (-ir+icos(Ai)J2„,+i(A)cos[(2m+l)X], 

F2rn = (-l)™sin(Ai)J2™(A)cos[2mX], 

G2™+i = (-l)™+icos(^i)J2™+i(^)sin([2m + l)X], 

G2ra = (-l)"sin(Ai)J2™(A)sin[2mX]. (B3) 



13 



Finally, we find for the coefficients of exp[Q{t) — Qi{t)] Here, Im is the modified Bessel function of the first kind 

of order m. 

Hm = I„i{Y) cos (mV) , 

Km = /™(r)sin(mT/), (B4) 

where we have introduced 
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